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A two dimensional hopping system with Rashba spin-orbit interaction is considered. Our main 
interest is concerned with the evolution of the spin degree of freedom of the electrons. We derive the 
rate equations governing the evolution of the charge density and spin polarization of this system in 
the Markovian limit in one-particle approximation. If only two-site hopping events are taken into 
account, the evolution of the charge density and of the spin polarization is found to be decoupled. 
A critical electric field is found, above which oscillations are superimposed on the temporal decay 
of the total polarization. A coupling between charge density and spin polarization occurs on the 
level of three-site hopping events. The coupling terms are identified as the anomalous Hall effect 
and the recently proposed spin Hall effect. Thus, an unpolarized charge current through a sheet of 
Ch ' finite width leads to a transversal spin accumulation in our model system. 

CN ■ I. INTRODUCTION 

The emerging field of spintronics tries to put the spin degree of freedom of electrons to use in a fashion akin 
J> . to electronics, where the charge of the electrons is utilizedJ^ One of the advantages of spintronics is, that spatial 
inhomogeneities of the spin distribution are not burdened with an energetic penalty, quite in contrast to inhomogeneous 
t , charge distributions. On the other hand, spin is not conserved, whereas charge is. This is a disadvantage of spintronics, 
] ' which means that spin polarization normally decays with time. Of special interest are techniques which influence the 
, spin by purely electrical means, without utilizing magnetic fields or magnetic materials. Spin-orbit interaction (SOI) 
provides such a mechanism. 

In recent years much has been done in this field for itinerant electrons, favoring very clean samples of high mobilityi^ 
. This investigation, on the other hand, studies the behavior of spins in a hopping system. Here, mobilities are very low, 
' either through disorder (Anderson localization), or through strong electron- phonon interaction (polaron formation). 
^ . Transport is thermally activated and incoherent. In the absence of effective spin scattering mechanisms, the low 
mobility might lead to long spin coherence times for a given spin coherence length, which is one motivation for this 
study, apart from the intrinsic scientific interest of the question of spin behavior in hopping systems. 

Previous studies of SOI in hopping systems centered on the influence of the spin dynamics on charge transport (e.g. 
Q ' magneto-conductance^iSiSiL2i2ii£), not on the spin dynamics itself. Furthermore, the effective spin coupling between 
• • hopping sites has normally been taken to be random(2i2iiiii2iiii4 which is appropriate for spin scattering on charged 
impurities, but not for SOI due to intrinsic or external electric fields which are constant over length scales of several 
hopping lengths. But such a mechanism would be required if one has the intent of affecting the spins in a systematic 
' way, i.e. realize spintronics. 

To be specific, we consider a system where the spin dynamics is solely determined by spin-orbit interaction through 
the Rashba-mechanisnii^iiSiiLiSi and no other spin scattering occurs (e.g. no magnetic impurities or spin scattering on 
charged impurities). The electronic system is two-dimensional (2D) and — in order to calculate explicit expressions 
— assumed to be ordered (small polaronic system). The disordered system is expected to behave qualitatively similar, 
but its investigation must be delegated to future publications. 

Two prominent examples of the interplay between charge and spin transport are the anomalous Hall effect and the 
spin Hall effect: The observation that a spin polarized charge current leads to a transversal Hall voltage, even without 
an external magnetic field, is called the anomalous Hall effect (see e.g. Ref. and references therein). The inverse 
effect, that, in materials showing the anomalous Hall effect, an unpolarized charge current leads to a transversal spin 
(but not charge!) current has been proposed by Hirsch^SiSiS and is called the spin Hall effect. Both of these effects 
will in the following be identified for our model system. 

This paper is organized in the following way: Sectionals concerned with the question of how to include the Rashba 
SOI into the hopping formalism. In Section IlIII the rate equations governing the evolution of the (1-particle) density 
matrix are derived. This is done by calculating the diagrams of second and third order (two-site hopping and three-site 
hopping) in the Konstantinov-Perel diagram technique in 1-particle approximation and applying the Markovian limit. 
The rate equations for the density matrix are transformed into rate equations for the particle density and the spin 
orientation, which is an equivalent representation, but much more lucid physically. It is found that two-site hopping 
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does not introduce a coupling between the equations for the particle (charge) density and the equations for the spin 
polarization. Three-site hopping processes introduce such a coupling and indeed lead to the anomalous Hall effect 
and the spin Hall effect for this hopping system. 

Sections |IVl and |V| give solutions of the rate equations for an ordered hopping system (i.e. a system of small 
polarons). First, bulk properties (i.e. no boundaries) are considered. A short account of the results of Sec. II VI has 
been published in Ref. |23- Then, in Sec.0 a strip of finite width is considered, which introduces boundary condition 
at the transversal edges. The stationary state with an unpolarized charge current shows spin accumulation at the 
boundaries, which is an expression of the spin-Hall effect. 

The appendices primarily deal with mathematical details. Appendix El considers the magnetic field due to spin 
accumulation as a possible means of detection of spin accumulation. 

II. RASHBA SPIN-ORBIT INTERACTION IN HOPPING SYSTEMS 

The first question which arises, is how to include the spin-orbit interaction (SOI) into the formalism of hopping 
transport. The Rashba Hamiltonian reads 

H^^^-a-{Kxp) + V{r), (1) 
Zm m 

where we introduce the quantity K = mae^/Ti^, which (i) has the dimension of an inverse length, (ii) is perpendicular 
to the two-dimensional plane (unit vector e^), and (iii) the length of which signifies the Rashba SOI strength a. The 
same expression Eq. ^ can be used for the generic SOI due to a spatially constant field E, in which case K = ^^^E. 
The symbol cr denotes the vector of Pauli spin matrices. 
Eq. can be transformed to 

H = ^{p-haxKf + V{r), (2) 
Zm 

while ignoring an irrelevant constant energy offset. Thus, in principle, the Rashba SOI can be dealt with as a SU(2) 
gauge potential. Since the gauge potential is a constant, this would be quite trivial, were it not for the fact that SU(2) 
is non-Abelian. If the gauge were Abelian, a wave function ^'(r) would be studded with the factor exp(icr • (K x r)) 
in order to include SOI. Due to the non-Abelian nature, there are additional higher order contributions in K. Those 
can be neglected, provided K is sufficiently small. 
The Hamiltonian which we consider reads 

mo- mva'aa' 

where H^^ph denotes the electron-phonon interaction and Hph is the Hamiltonian of the phonon system. Here, aj^^. 
(omcr) is the creation (annihilation) operator of an electron at site m in spin state cr, is the site energy, and J^j/*^ is 
the resonance integral between the states ma and m'a' . Without SOI, the resonance integral is diagonal in spin space 
and will be denoted as Jm'm in the following. Note, that the usual spin splitting due to Rashba SOI is not manifest 
here, since Kramer's degeneracy demands, that the two spin states of a localized eigenstate have to be degenerate. 

The central question for the inclusion of SOI into the hopping Hamiltonian now consists in determining the spin 
structure of J^/*^. The four numbers belonging to each pair of sites m'm are collected into a 2 x 2-matrix in spin 
space Jm'm- Clearly Jm'm ~^ Jm'ml2 for K 0, where I2 denotes the two-dimensional unit matrix. Since the 
Hamiltonian has furthermore to be Hcrmitic and invariant under time-reversal symmetry, the spin structure is given 
by an SU(2) niatrix^Si^iiiiiLSiiSSiSSiSi Jm'm — exp(— icr • Am'm)Jm'm, where all components of the vector A have to be 
real numbers and the condition A^'m = —-^mm' must be obeyed. The arguments after Eq. Q show that the term 
linear in K reads —icr • (K x Hm'm), where Rm'm = Rm' — Rm is the distance vector between the sites, and Rm is 
the coordinate vector of site m. Thus, 

'-'m'm — ^ '-'m'm^ V 7 

For this expression for J to be valid, K has to be sufficiently small, as explained in the paragraph after Eq. In this 
case, this means, that the condition K x Rm'm ^ 1 must be valid for the distance vectors Rm'm, which are relevant 
to the case at hand (e.g. the typical hopping length or the lattice constant). The same expression, but applied to the 
Green's function of a hopping system, has been derived in Ref. ^ 
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In some sense, this is equivalent to the Holstein transformation's, -where the influence of a magnetic field on a hopping 
system is taken into account through a phase factor, and where higher order (in the magnetic field) corrections are 
neglected. 

We are now prepared to derive rate equations governing the hopping system with the inclusion of Rashba SOI in 
the next section. 



III. RATE EQUATIONS 

We apply the Konstantinov-Perel diagram technique in order to obtain rate equations for the density matrix. The 
off-diagonal elements of the (1-electron) density matrix {al^ic'Cima) can usually be neglected in a hopping system 
(otherwise, the transport mechanism would not be hopping). Here, we have to keep in mind, that the correlations 
between different spin states on the same site contain information regarding the (expectation value of the) spin 
orientation. Furthermore, the model which we consider, aims at systems where the spin degrees of freedom do not 
de-cohere (lose their phase memory) while the electron stays on a site. Thus, we must retain the off-diagonal elements 
in spin space on the same site (i.e. diagonal in the site index) of the density matrix. In this way, by including the 
SOI the occupation probability pm = {'An'^rn) becomes a 2 x 2-matrix Pm\a'a = {o^ma' '^ma) m spin space. 

In order to adapt the Konstantinov-Perel diagram technique to the SOI case, each site index can be thought of as 
additionally containing the spin index. The summations run over all sites and also the two spin states. The electron- 
phonon interaction is unaffected by the spin state. Thus, the extension of the formalism reduces to (i) studding each 
site index with a spin index (and extending the corresponding sums) and (ii) allowing expectation values which are 
non-diagonal in spin space. The number of terms in a given diagram proliferates quite rapidly with the order of 
the diagram and makes its computation tedious. The calculations are greatly simplified, if one replaces the explicit 
summation over the spin indices by matrix multiplication in spin space. This can be achieved by collecting the 
interaction matrix elements J (which depend on two site and two spin indices) into appropriate matrices in spin space 
Jm'm (which only depend on two site indices) and doing likewise with the density matrix elements. 

The second order diagrams (two-site hopping) yield in one-particle approximation and in the Markovian limit the 
rate equations 

— o _ V le"''^ (^''^'"i'")o i<T (KxR„.i,„)w - n W \ fS^ 
mi 

where the transition rates Wm'm are the same as are obtained without SOI. 

The physical meaning of Eq. ISJ is easier to assess, if the following transformation is applied: Representing the 
2 X 2-density matrix at site m by pm = \pmI2 + \pm ' where I2 is the two-dimensional unit matrix, the occu- 
pation probability = Tr(pm) and the (expectation value of the) spin orientation = Tr(<T/5m) are introduced. 
Please note, that in the following we will use the wording "spin orientation" without emphasizing each time, that its 
expectation value is meant. Of course, a spin-1/2 particle does not have a classical spin orientation. 

The factors (matrices in spin space) surrounding pmi in Eq. Q are inverses of each other. Thus, by taking the 
trace (in spin space) of Eq. @ these factors disappear and one obtains 

^Pm = ^ {PmiWm^m ~ PmWmnn} ■ (6) 
mi 

On the other hand, multiplying Eq. Q by cr and then taking the trace, one obtains the rate equation for the spin 
orientation (see Appendix fX)) 

^Pm ^ ^ "^-D^jifj • PjYii^mim Pm^mmi^ j C^) 

771 1 

where the quantity l) is a 3 x 3-matrix describing a rotation about the axis Amim = K x R„jj„j. The Eqs. ^ and 
{Tj) are the sought for rate equations which replace Eq. ||SJ|. 

One can see, that in two-site approximation the equations for the occupation number and the spin orientation 
decouple. Note, that the terms of order are retained in Eq. iQ, in addition to the terms of order K, since the 
third-order terms, to be derived next, give contributions only to second order. The behavior of an ordered (polaronic) 
hopping system obeying Eqs. © and Q is considered in Section Hvl 

Next, the three-site probabilities (third order diagrams) are calculated. These contributions yield higher order 
corrections to the two-site probabilities, as well as some terms, which go beyond the physics of the two-site expressions. 
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Only the latter type of third order terms will be retained in the following calculations. The products of spin matrices 
occurring in the third order terms are calculated in Appendix ^ The rate equations up to terms of order and 
neglecting superfluous terms as mentioned above read 



dt 



} (8) 



mim2 

}• (9) 

mim2 

Eqs. © and © are the generic hopping rate equations for the model considered here (only Rashba SOI affects 
spin). The transition rates Wm'm, W^mi?n2m; ^^'^ ^rnxm^m ^I'S independent of spin and are determined by the type of 
hopping considered (e.g. small polaron hopping or hopping between impurities). In the remaining part of the paper, 
we focus on an ordered (polaronic) system, but Eqs. © and (PJ would also be the starting point for an investigation 
of a disordered system. Summing Eq. (jS)) over all sites m yields ^ ^.^^ pm = 0, thus, the derived rate equation for 
the particle density obeys particle (and charge) conservation. 

Whereas the second order transition rates Wmim are those of a system without SOI, the third order transition 
rates W^^^^^ differ from the corresponding quantities without SOI. Without SOI, the third order transition rates 
(let us call them W^^'^ summarily) are important for the Hall effect. The W'^^^ are divided into parts which are 
symmetric in the magnetic field and parts which are anti-symmetric. Only the latter are important for the Hall 
effect, thus the symmetric parts are normally neglected. Here, since the magnetic field is zero, the anti-symmetric 
part vanishes. Thus, we have to calculate the usually neglected symmetric parts of VF'^^. Furthermore, both, real 
[W^) and imaginary iW^) part of W^'^^ find entry in the formalism, whereas only the real part of the second order 
W occurs. Nevertheless, is connected with the Hall mobility of a hopping system. The calculations are shown in 
Appendix 151 

If the system is ordered (i.e. polaronic), it is advantageous to transform the equations into wave- vector space. 
Note that, since the spatial vectors are two-dimensional, the wave-vectors q will also be two-dimensional. This 
transformation leads to the rate equations 

|p(q) = [l^(q) - Wm p(q) + [W^(q) - W«(0)] • p(q) (10) 

and 

^p(q) - (W{ci) - l^(0))p(q) - 2#(q) • p(q) - 2(K x W(q)) x p(q) 
at 

+ (W«(q) - W^(0))p(q) - (W^(q) - W^(0)) x p(q) (11) 

Here, 

W(q) - Y^W^(q) (12) 



(13) 



W«^(q)=KoK. ( Ax A)M^i^./(qq^) 



(14) 



qi=0 



where the last equation is to be understood as applying to the superscripts R and /, respectively. 
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The basic transition rates are 



iy(q) =^e"i-^'"M^Om (15) 



M^'''(qqi) = E e^^'^'-e^i-^-i^^^^^ (16) 

mim 

In an isotropic system, taking the long wavelength limit, they read 

Wiq) = WiO) - Dq^ - i^q • E (17) 

and 

W^^(q) =iKoq. (E X K)W^^ (18) 

the latter being calculated for a triangular lattice (see Appendix Ib)| . Here, D is the diffusion constant, /i the mobility, 
and the quantities and are constants defined through Eq. (|B8|I . 
The rate equations thus take the shape 

|p(q) = -iq-j(q) (19) 
with the particle current density (proportional to the charge current density ej) 

j (q) = MEp(q) - iDqp{q) + W'^K X E o K • p(q) (20) 

and 

^p(q) = -iq • J(q) + Q(q) (21) 

with the (tensorial) spin current density 

J(q) = ^E o p(q) - iZJq o p(q) + 4i:i/3K • p(q) - 4i:'p(q) o K 

+W'^K X E o Kp(q) - W'K x E o K x p(q) (22) 

and the spin source density 

Q(q) = -4DKV(q) - 4DK o K • p(q) + 2^(K x E) x p(q) (23) 

Eq. H19|l is the continuity equation for the particle (charge) density. Eq. (|21|l contains a source term in addition 
to the divergence of the current, since there is no conservation law for spin. The first two terms contributing to the 
charge current H2U|I are the drift term in the electrical field and the diffusion. The third term only occurs, when there 
is a finite z-component of the polarization (i.e. a magnetization M^). The resulting current is directed perpendicular 
to the electric field. Thus, this term corresponds to the anomalous Hall effect. The corresponding mobility is seen to 
be iJ,yx = W^K"^. The quantity /i^x/M is related to the Hall mobility (see Appendix Ib)| . 

The spin current Eq. (|22|) mainly consists of the drift term and some diffusion terms. But the most interesting 
contribution is the fifth term on the right-hand side, since it is also present for vanishing spin polarization p — 0. 
This term describes a current of z-spins into the direction perpendicular to the electric field (K x E). Thus, it is an 
expression of the spin Hall effect. The co-efficient is the same, as the one responsible for the anomalous Hall 
effect. 

In order to simplify notation, dimensionless quantities are introduced: dimensionless space x — Kr, wave vector 
^ = q/K, time r = DK^t and electric field e = fiE/{DK). The transformed third order probabilities are denoted 
as eji = ^y-W^ and e/ = ^^W^ . Note, that en can be expressed as ^yxc/ ^. Thus, cr is an analog of the Hall 
mobility multiplied by the electrical field. Furthermore, the co-ordinate system is fixed such, that oc K, e^; oc E, 
and Gj, cx K X E. Then K = Ke^ and E = Ee^. This transforms the rate equations H19|) and (|21|l to 

r) ^~{e+ i^.e)p(l, t) - ieR£,ypM. r) (24) 



6 



and 

t) = -(1^ + iCx£)p - 4i(e^ X ^) X p - ien^yC^p + ie/C^e^ x p 

-4p - 46^^^ + 2eey x p (25) 

The three terms on the right-hand side of Eq. (|24|) are the diffusion contribution, the drift in the electric field and 
the anomalous Hall effect. The first terms on the right-hand side of Eq. (|25|l are: diffusion, drift (Ohmic current), a 
precession term for diffusing spins, the spin Hall effect, and a precession of the spin about the z-axis. The last three 
terms contribute to the spin decay. The first two of those can be expanded as —Ae^Px — '^^yPy ~ Sie^pz, thus the decay 
rate for the z-component is twice as large compared with the in-plane components. 



IV. BULK SYSTEM 



By setting the wave vector ^ to zero, the evolution equations for the quantities integrated over the whole system 
(total charge and total spin polarization, e.g. p{^ = 0, r) = J (Pxp{:x., r)) are obtained 



— p(0,r) = 



-4p- 4e2p^ -I- 2eey x p. 



(26) 
(27) 



Eq. H26|) is an expression of charge conservation, ft is to be seen, that the evolution of the space integrated quantities 
does not depend on the third-order parameters e^j and e/. 
The three components of Eq. H27|l read 



d 
d 



'Ap^ + 2£pz 
-8pz - 2epx, 



the solution of which is 
Px{t) 



pM cosh(2ryT^) + /^^W + ^^^W sinh(2rv/T^) 



Py{r) = e-'^Py{0) 



pAr) 



-6t 



pM cosh(2TyT^) - sinh(2ryT^) 



(28a) 
(28b) 
(28c) 

(29a) 
(29b) 
(29c) 



Px(0), Py(0), and ^2(0) are the initial conditions. Note, that there is a critical electrical field strength €c = 1. In 
smaller fields the total spin polarization decays exponentially with time, but in larger fields, it additionally rotates 
with time. This critical electric field Ec = DK/ p, can also be expressed as Ec — kBTK/e, provided the Einstein 
relation between D and p is valid. Furthermore, without electric field, the z-component of the total spin polarization 
decays twice as fast as the other two componentsi^iSS 

Taking for definiteness the initial condition p(0) = , one obtains 



Px{t) = 
P,(t) = 



-6t 



: sinh(2Tv/l - £2) 



-6t 



COS. 



h(2TVr 



sinh(2T\/f - e2) 



(30a) 
(30b) 
(30c) 



for e < 1, whereas for e > 1 the hyperbolic functions become trigonometric functions. Expressed differently, for e > 1 
the eigen- values of the set of equations 128|l become complex, which means that its eigen-modes become oscillatory. 
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Thus, 



P,(t) = e-^^^^sm(2Tv/e2-l) (31a) 
P,(t) - (31b) 



P,(t) = e''" __ cos(2tV£^-1 + 0), (31c) 
ve — 1 



with sin0 = 1/e. Note, that the frequency of the polarization oscillation 2%/? — 1 increases with the electric field 
e. Apart from the decay factor e~^^, the polarization plr) describes a slanted ellipse in the x-z-plane, i.e. the plane 
spanned by the Rashba field K and the electric field E. The co-ordinates x — p^e^'^ \/ ~ ^'^^ ^ ^ Pz^'^^^ ~ 1/e 
determine the ellipse + 2xz sin + = cos^ 0. The inclination of this ellipse (angle between z-axis and the main 
diagonal) increases with electric field. For e 1 the angle — > 7r/2, i.e. the ellipse degenerates to a line with an 
inclination of 45° to the perpendicular direction. In the limit of strong electric field e — > oo, the angle vanishes, and 
the ellipse becomes a circle. In order to observe these oscillations, the corresponding frequency should be larger or of 
the order of the inverse decay time constant, which leads in this case to the approximate condition e ^ 10. 

Let us consider next the case of an inhomogeneous initial condition. We take a z-spin starting at the origin 
p(x, r = 0) = (5(x)ez. The total spin polarization is still described by Eqs. 129|l . In order to obtain the local spin 
orientation, one has to solve Eqs. 1)24(1 and H25|l . This is done by solving these equations in Laplace-space and the 
following back-transformation from Laplace-space s to time r and wave-vector space ^ to x. 

The asymptotic behavior of the local spin polarization for large times r — > oo and small distances |x| ^ l^fr (Large 
distances will not be relevant here, since they are strongly suppressed by the "diffusion factor" exp(— x^/4t).) is 

r (x + er)2 7 1 3 , /VTS, A 
P.(x,r)=exp| r|^Jo(^ — (32) 

Note, that in this approximation (long times and short distances), the sign of the polarization — determined by 
the argument of the Bessel function Jq — only depends on the distance to the origin, and is furthermore not time 
dependent. This spatial oscillation is due to the precession terms in the rate equations. The exponential factor 
exhibits a decay term with (dimensionless) time constant 4/7 and a drift-diffusion term with drift velocity — e and 
diffusion constant 1. 



V. STRIP OF FINITE WIDTH 



The distinguishing feature of a strip of finite width is, that there is no current (neither charge, nor spin) across 
the transversal boundaries. Thus, the spin-Hall current must be compensated by a diffusion current, i.e. a finite 
spin density can be expected to occur near the boundary. Indeed, the following will show, that a transversal spin 
accumulation will occur in this situation. For a broad strip, the non-zero spin density is mainly constricted to near 
the boundaries. 

The strip is oriented, so that the x-direction is along the length of the strip (the direction of the current flow), the 
y-direction is along the width of the strip and the z-dircction is perpendicular to the 2D strip (see Fig. The strip 
extends from y — —6 to y = b, thus, its width is 26. Going from wave-vector space back to x-space and assuming that 
p(x) and p(x) = (p^, Py, Pz) depend only on the y-coordinate, the following set of differential equations determine the 
stationary state 

p" + eRp', = (33a) 
p'^ + eip'y + 2€p,-Ap, = (33b) 
p';^-eip[^+Ap',-Apy = (33c) 
p'^ + eRp' ^Ap'y-2€p,-Sp, = (33d) 

where the prime denotes the derivation with respect to y. The boundary conditions at y = ±6 are obtained by 
equating the normal component of the currents H2(J|) and (|22|1 at the boundary to zero. This yields 

{p + eRPz)\v=±b = (34a) 

{p'x + ^iPy)\v=±b = (34b) 

{p'y-tiPx + '^Pz)\v=±b = (34c) 

{p', + eRp - Apy)\y=±b = 0. (34d) 



8 



One has to keep in mind, that we want the solution to lowest order in K. The zero-th order solution of Eqs. H33a() 
to (|34dp is obtained by removing all terms at least linear in K, i.e. setting ej — and e^j — 0. This immediately yields 
the solution p(jj) = const, and p{y) = 0. This solution is inserted into the terms containing en and e/, which yields 
the equations for the first order solution 



p" = 

p'^ + 2ep,~ip^ = 

p'^ + Ap',-Apy = 

Ap'y - 2ep^ - 8p, = 



(35a) 
(35b) 
(35c) 
(35d) 



with the boundary conditions 



ip%=±b = 

{Px)\y=±b = 

{Py + 4.p,)\y = ±b = 

(p'z ~ 4:Py)\y=±b = <^RP- 



(36a) 
(36b) 
(36c) 
(36d) 



The solution of this set of equations can easily be constructed, since the only difference to the equations of zero- 
th order is the finite right-hand side in the boundary condition l|36d|l . Again, p{y) — const., and the 3 vectorial 
p-components are weighted sums of functions sinh(A?/) and cosh(Ay) where A takes on the three values: Ai = 2, 
A2 = V-2 + 2iVe2 + 7 = VV8 + e'^ -T + i VVsT? + 1 and A3 = a/-2 - 2i\/"e2 + 7 ^ x*. Thus, to lowest order in 
K, 



Pxjy) 

<^RP 

Pyjy) 

^RP 

p^y) 

£rP 



16e 

- — A2C2S3 sinh(2y) + — ( I6C1S3 -|- — A3S1C3 ) sinh(A22/) + 



4e^ 



A2C2S3 cosh(2j/) 



Ai-4 



I6C1S3 



2e 

a 

4A2 

a 



c.c. 



I6C1S3 -I- — A3S1C3 ) cosh(A2?/) -I- c.c. 



A3S1C3 sinh(A2y) + c.c. 



(37a) 
(37b) 
(37c) 



with the abbreviations 



A2C2 



A2PA2S1C3 + 16{Xl + 12 + e^)ciS3 



(38) 



and 



Cj = cosh(Ai6), Si = smh{Xib) (i = 1,2,3). (39) 

One can see, that the electron density p and the spin Hall co-efficient only influence the amplitude of the po- 
larization. They do not affect the spatial dependence of p{y). The x- and z-component of the polarization are 
anti-symmetric with respect to a sign change of y, whereas the y-component is symmetric. 
For a large electric field e ^ 1 

pAy) « -sgn(y)^ sin (Vi(6 - |y|) + J) e--^(^-l^l) (40a) 

Pyiv) « -^sin(Vi(6-|z/|))e-^(^-l^l) (40b) 

p.iy) « -sgn(2/)^ cos (^^(6 - |y|) + J) e'v^^^-I^D (40c) 

The main effects of increasing the electric field e consists in a sharper concentration of the polarization near the 
boundary (due to the exponential term) and an increase of the amplitude of the x- and z-component of the polarization 
(because Cfl/Ve = V^P-yx/p- oc ^/e). 

For 6 3> 1 (i.e. in natural units: width 3> K~^) the spin polarization is significant only near the boundaries due to 
the hyperbolic functions. This suggests that there is a "spin diffusion length" governing the spatial extent of non-zero 
p near the borders. This suggestion can be fleshed out in the following way: We consider the half plane y >0, set the 



9 



electric field e to zero and apply the boundary condition p{y = 0) = Po- The differential equations for the stationary 
state (and homogeneous in x) are 

pL'-4p, =0, p'; + Ap',-Apy = 0, p^' - - = 0, (41) 

which conceptually corresponds to an equation with the structure ZJAp — p/t = 0, only that here the decay rate 
1/r depends on the direction (thus, it is a tensor), and there is an additional "torque" term. Since in Eq. H41() all 
co-efRcients are numbers (in dimensionless units) the resulting "spin diffusion length" ^/Dt is also a number of the 

order of one (in the present case: 1/2 for the x-component and 1/Vv8 — 1 ~ 0.74 for the other two components, 
corresponding to the real parts of Ai, A2 and A3 at e = 0). Thus, expressed in natural units, this length is times 
a constant of order one, i.e. this spin diffusion does not introduce a new length scale in addition to the Rashba length 
scale K^^. 

The spin current of a non-polarized electrical current in an infinitely extended plane consists of 2;-spins being 
transported perpendicularly to the electrical field E (see Eq. (|22|) '). Here, this spin current is "transformed" into spin 
accumulation by applying boundary conditions. Seen in this light, it is perhaps somewhat surprising that the resulting 
spin polarization has finite x- and y-components, in addition to the z-component. But this can be explained by the 
following observation: The calculated spin polarization p{y) is — though stationary — not static. It represents an 
equilibrium between the "spin supply" by the spin Hall effect on the one hand, and other contributions of the spin 
current as well as the spin decay, on the other hand. Even though the "supply" only yields z-spins, the tensorial 
nature of the current leads to spin precession, such that in the equilibrium state all components occur. 

Due to the spin polarization (magnetization) a magnetic field H7\/ is generated, which may in principle allow to 
detect this polarization. In order to actually measure this magnetic field, one has to take into account, that there is 
also a magnetic field H/ due to the electric current through the strip. These fields are calculated in Appendix lO 

In principle, the magnetic field Hm + H/ interacts with the spin degrees of freedom through the Zeeman energy 
and also applies a torque to the spins. These effects have been neglected in this paper. 

The (in-)significance of the back-reaction of the magnetic field due to the spins on the spins can be estimated in 
the following way: Using the magnetic susceptibility x the Bohr magneton /is, the spin polarization due to the 
back reaction can be expressed as p^^ — xH/fJ-B ~ X^/ f^B = XP- Thus, the influence of the magnetic field on the 
spins can be neglected, provided the condition x ^ 1 is fulfilled. 

VI. DISCUSSION 

We have studied the spin dynamics of a hopping model where Rashba SOI is the sole spin scattering mechanism. 
Rate equations have been derived in the one-particle approximation and the Markovian limit, including two-site and 
three-site hopping processes. Further calculations are performed for the case of polaron hopping, where one can work 
in wave- vector space. 

As expected, an initial spin polarization has been found to decay exponentially with time, if the electric field is low. 
But, there is a critical electric field, above which the total spin polarization shows oscillations (rotations) overlaying 
the decay. This can be understood as follows: 

The rotation of a spin initially at a certain site mi during a single hop to another site m2 depends only on the 
difference vector Rm^mj . The spin rotates in the plane spanned by the vectors K and Rmami • Let us consider a spin 
starting at the origin and being orientated along the z-direction at i = 0. The motion of the electron in the radial 
directions leads to concentric "waves" around the origin of spins having the same z- and radial component. But, due to 
the diffusive motion, the electron does not only take the straight path radially outwards, but explores an increasingly 
entangled trajectory, thereby randomizing the spin orientation. The random motion only affects a decrease of the 
magnitude of the polarization (it neither creates polarization, nor does it favor any specific orientation). Thus, the 
preferred orientation given by the radial motion is preserved, but its magnitude diminished. 

An electric field does affect the charge transport, but in two-site hopping approximation it does not affect the 
elementary processes responsible for the transport of spin orientation. Thus, in this case, the spin polarization is 
simply the product of particle density and spin orientation, which has a constant magnitude itself (quite in contrast 
to the situation with three-site hopping processes). In this way, one can see, that the electric field does not change the 
spatial distribution of the spin orientation (e.g. the sign of the z-component), but only changes the temporal evolution 
of the spatial charge distribution (the "weight" of the local spin orientation) . This explains the time independence of 
the argument of the Bessel function in Eq. (|32|l . 

On the other hand, in a very large electric field an electron follows preferably a straight trajectory along the field 
(large drift component in comparison to the diffusion). The importance of diffusion about this trajectory is reduced, 
thereby also reducing the randomization of the spin. This explains, why it is possible for the electric field to become 
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so strong, that the spatial oscillations of the spin orientation survive the averaging over the whole volume, so that 
they appear as temporal oscillations. 

Without an electric field, the decay of an initial spin orientation is found to behave anisotropic: The z-component 
(perpendicular to the plane) decays twice as fast as the in-plane components. 

The inclusion of three-site hopping terms leads to the occurrence of a transport co-efficient (e^ or fiyx) which provides 
a coupling between charge and spin transport. A spin current produces a transversal charge current (anomalous Hall 
effect), and a charge current produces a transversal spin current (spin Hall effect). In a two-dimensional strip of finite 
width, this transversal spin current manifests itself as a spin accumulation at the boundaries of the strip. 

Finally, some remarks about the experimental relevance are in order here. The Rashba SOI strength is widely 
given as being of the order a « 10~^ eV cmi^ti^ It is possible to change this value by applying a gate voltage to 
the two-dimensional planed, where a variation by at least 50 % has been shown to be realizable. This might be 
conceptually favorable for the implementation of spintronics devices, where electric means to control spins are sought. 

Using the value of a mentioned above, one obtains K « 1/(76 A), assuming that the effective electron mass is 
equal to its elementary mass. Thus, the condition K x Rm'm <C 1 for the validity of our theory demands in this 
case, that the typical hopping length is shorter than 76 A. This value is conveniently large, so that the range of 
polaronic nearest neighbor hopping is well within the validity range. But impurity hopping or variable range hopping 
are also not out of the question. Furthermore, for a smaller Rashba co-efRcient, the permitted length increases and 
it is furthermore modified by the effective electron mass. If the Einstein relation between diffusion co-efficient and 
mobility is assumed to be valid, the critical electric field Ec = ksTK/e for the oscillatory behavior of the total spin 
polarization is estimated to be about 10^ V/cm at a temperature of 100 K. 

The calculations in Appendix^ show a close connection between Hall mobility and the co-efficient for the spin Hall 
effect. Thus, a large Hall mobility favors also the occurrence of the spin Hall effect. We take Ref. |^ as an example. 
There, the conductivity and the Hall mobility of the hexaboride compounds Eui-ajCaa^Bg are measured, and the 
authors conclude, that the transport mechanism is polaron hopping (in a certain range of concentrations x and for 
high temperatures). The reported Hall mobility in the polaronic hopping regime is of the order of 100 cm^/(Vs). 
Using this value, Eq. ljB12|l gives the estimate fiyx/lJ- ~ 0.23, i.e. the transversal mobility for the spin Hall effect is 
about a fourth of the longitudinal mobility. Taking Eq. I|37c|l into account, this gives a relative spin polarization 
(z-component) at the boundary of a two-dimensional strip \pz(b)/p\ of about 5% at the critical electric field strength 
e = 1. 
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APPENDIX A: CALCULATION OF THE SPIN MATRICES 



Applying the identity (cr • A)(cr • B) = (A • B) -I- icr • (A x B), valid for any vectors A and B, to the series expansion 
of the exponential function, one easily obtains the relation 

e'-^-^^coslAI+ia-A^^, (Al) 

|A| 

which is valid for any vector A. 

The quantity Tr (cre^"^-^™i'"/5mj^e"^^'"i'") occurs in the calculation of the contribution of two site diagrams to the 
rate equations. It can be transformed to Tr spin factors in 

the trace can be calculated to be 

g-A„,„^g-i^A„,„ ^ D^^^a, (A2) 

where the quantity D is a 3 x 3-matrix describing a rotation about the axis Amim = K x Rm^m and through the 
angle 2|Amj,„|. This is valid to quadratic order in A. Explicitly, 

Dmirn = 6^ X o X + siu |2A|(e/i o - o e/j) + cos |2A|(e^ oe^ + eRO bb)- (A3) 

The symbol "o" denotes the dyadic product and the orthogonal unit vectors are = Ji/K and e^; = Rmim/IR-miml- 
Here, it is used that K is perpendicular to the plane, i.e. K • R = for all vectors R within the plane. Retaining in 
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this expression only terms of up to second order in K (since ^ ^ 1), one obtains 

Drmm. ^ I3 + 2|A|(eK o - o e^,) - 2\A\'^{e^ oe^ + eRO br), (A4) 

where I3 is the identity matrix in three dimensions. Note that det{D) = 1 + 0{A^), where 0{x) is Landau's big-0 
notation, i.e. that the length of the vector p does not change (up to 0{A'^)) due to the application of the operator 
D. Omitting the A^-term, the length of the vector would change in second order of A; then, D would not represent 
a rotation. 

The products of spin matrices in the three site diagrams are calculated using the following scheme: After rotating 
the operators within the trace, so that p is the last operator, the spin factors preceding p in the trace are calculated 
in order to obtain an expression of the form / + gcr. Here, / and g are scalars in spin space. Thus, taking the trace 
yields Tr ((/ + gcr)p) = fp + gp, an expression where the spin operators have disappeared and are replaced by the 
particle density and the spin orientation. Using the operator identity e^e^ « qA+b^ia,b]/2^ where the error is of third 
and higher order in the operators A and B, one obtains 

and 

g-io-A„„2g-i<T-A„2„.jg-i<T-A„j„ gio-.(A„„jXA„„2)_ ^^g-j 

Using Eqs. ljA2|l . ljA5|l . and (jASp . which are all valid to second order in A (i.e. also in K), and Eq. ljAl|l . the spin 
structure of the three-site terms can be calculated, e.g. 

Tr(e---A„™,g-i-A,„,„,g-i^.A,„,„^) « cos|A2|p + i^^^^j^A2-p 

IA2I 

« p + iA2-p, (A7) 



sin|A2| sin|A2| 
cos IA2IP + 1 ' I A2P ' I A2 X p 



~ i^mam ' P + - Aj X p, (A8) 

where the abbreviation A2 ~ Am„ii x A„„„2 is used, and the second approximation for each relation is the truncation 
to second order in K. 

The expressions (|A7|) and (|A8|) occur together with the (complex) third order transition rate W^^\ Adding the 
corresponding complex conjugate diagram, one obtains some terms which contain the real part = 2Ile{W^^^) and 
some which contain the imaginary part = 2Im(W^'^^^). Without SOI, only the real part is relevant, since there is 
no complex pre-factor in this case. 

APPENDIX B: THIRD ORDER TRANSITION PROBABILITY FOR SMALL POLARONS 

The third order diagrams (three-site hopping) give the additional contribution to the right-hand side of the rate 
equation (O 



i ^ I - e-'- A™™2 e-i- A^^-M p,„^ e-'- A'-i" 



mim2 



4_ „-i<T-A,„„, ' -i<T-A„,„i, -iCT-A„2„Tj/(3)* 

- Ame-'*"-^'""! e-*- A-i™^ e-i- A™„„ 

^g-i^.A„„,g-i^.A„„„,g-i^.A„,,„^^^yf^052^^^^^J _ (Bl) 



The third order transition rates are Wj^^^^„^ = wL!m2m + W4iim2m = 2Re(WL!m2rn) and W^^^^^^ = 

(3) fS)* (3) 

-i{Wmim2m - Wmim2m) = 2Im( Wj^/m^m ) , derived from the quantity 

^mlm2m ~ Jmm2Jm2miJmimI , (B2) 
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where the symbol / denotes the integral 

/•oo 

/(E) = e"^'^^ / dti I dt2 I exp 



OO 







2sinh(ntj,/3/2) 



/3 /? (3 

COSUJq{tl + 1^-2) + C0SU)q{t2 + 1^2') + C0SW5(il — t2 + ^^'2^ 

X exp j ^eE • R„,„iii - -^eE • R™„2*2 ^ (B3) 



with the abbreviations 



and 



^ l7(g)P(l-cosq.g) 
q 

The sums runs over the phonon wave- vectors q. 7(9) is the electron-phonon interaction constant, LUq the phonon 
frequency, N the number of phonon states, (3 = l/(fcsr), and g is a vector connecting two nearest-neighbor sites. 
Note, that the integral is independent of the direction of g. Thus, the resulting expression is isotropic. 

Taking two times the real part, the integral over t2 would extend from —00 to 00. This would be the expression to 
be calculated for the small polaron Hall effecti^^i^ Here, we cannot take this step, because we also need to calculate 
the imaginary part of the expression ljB3|l . 

The three-site rate obeys the principle of detailed balance 

^(3) =^(3) -/3eE.R„„, /ggN 

and has the further symmetry 

Using these symmetries, the long wave-length limit of the third order term Eq. H14|l for a triangular lattice (hexagonal 
crystal symmetry) can be calculated as 

W«(q) = iK o q . (E X K)|^Re(<^^). (B8) 

In order to calculate the imaginary part of W^'^^ has to be used instead of the real part. The indices 0, 1 and 
2 denote the sites located at the vertices of an elementary triangle of the lattice. Wqi2 has to be calculated in zero 
electric field. Thus, there only remains to determine the integral /(□) in order to obtain the third order coefhcients. 

For a sufhciently large phonon dispersion, the main contribution to the integral / comes from a saddle point (^i^, i2s) 
at (i2?i/3/3, i?i/3/2). Under the condition hujq ^ fc^T (high temperatures compared with the Debye temperature), the 
exponent in the integrand can be expanded in a series in ti and <2, which is truncated after the second order terms. 
The integral can then be calculated and one obtains 



^"bt 1-iW — —e-'^BT ^ (B9) 



where 



Ea-T.^- (BIO) 
q 

Combining Eq. I)B9|I with Eq. (|B8p . the mobility fiyx connected with the anomalous Hall effect and the spin Hall 
effect reads 
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where no = eo? /Ti and a is the distance between nearest neighbors (lattice constant). This expression can be compared 
with the Hall mobility fiH for small polaron hoppingSSr to give the relation 

^ = 2(Xaf^ = ^^^. (B12) 

The effective activation energy for the processes (e.g. the spin Hall effect) is 4£'a/3 [see the real part of Eq. 
(IB9|l ]. which is identical to the situation for the ordinary Hall effect. On the other hand, the activation energy for the 
processes is Ea, i.e. ^ for small temperatures fesT <C Ea- 



APPENDIX C: MAGNETIC FIELDS DUE TO SPIN ACCUMULATION AND ELECTRIC CURRENT 

The magnetic field generated by the spin polarization (magnetization) M(r) = M(rj,) = fispiy) in the strip reads 

27r rJ + [y- y'y 

Note, that the field Ha/ has no z-component, even though the spin polarization has a finite x-component M^. The 
magnetic field Hm may be exploited in order to detect the spin accumulation. In this case, one has to keep in mind, 
that the electrical current in x-direction also leads to a magnetic field H/.-^° Both fields have different symmetry. 
Thus, if Hm is not very much smaller than H/, a measurement of the magnetic field near the strip allows a detection 
of the spin accumulation. 

Explicitly, the magnetic field generated by the (homogeneous) current / = eb^Ep within the strip reads 

H/(y, z) = -rd arctan , + ^ log -^^^^ ■ (C2) 



Anb \y^ + z'^-b^J 8nb "^yz^ + iy + b)^ 
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FIG. 1: The co-ordinate system chosen for the strip of finite width. The charge current within the strip is non-polarized. But 
at the transversal boundaries y = ±6, a finite spin density occurs due to €r, i.e. this situation leads to spin accumulation. 
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FIG. 2: The magnetic field due to the spin accumulation in the 0-?/-plane. The parameters are 6 = 1, e = 1. The horizontal 
line denotes the extent of the strip. 
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FIG. 3: The magnetic field due to the current through the strip. The parameters are the same as in Fig. |5| 



